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Abstract. IRAS 13349+2438 was the first quasar suggested to 
host a warm absorber with internal dust. We obtained new HRI 
data for this object and derive constraints on the X-ray emitting 
region of IRAS 13349+2438 by investigating its X-ray extent 
and lightcurve. Rapid X-ray variability is detected with a factor 
of two change in count rate within one day and 20% variability 
within about three hours, which practically rules out any signif- 
icant contribution of scattered X-rays. We present for the first 
time a detailed, self-consistent modeling of the ROSAT PSPC 
and ASCA spectrum of IRAS 13349+2438 in terms of a warm 
absorber including dust. It turns out that both a dust-free and 
a dusty single-component warm absorber fail to explain the 
ROSAT, ASCA, and optical data simultaneously. We discuss 
possible explanations such as a variable and/or complex warm 
absorber, a variable soft X-ray continuum and instrumental ef- 
fects. 
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1. Introduction 

IRAS 13349+2438 is the archetypal infrared-selected quasar, 
discovered in a sam ple of bright IRAS 12/x sources by Beich- 
mann et al. ( 1986| ). It is optically bright (my = 14.7) and 
located at a redshift of z = 0.10764 ± 0.00027 (Kim et al. 



1995 ). Broadband investigations revealed that the spectral en- 
ergy distribution is dominated by infrared emission. Its bolo- 

erg 



metric luminosityf]bet ween 0.34 and 100/i is w 2.3 x 10 46 • 



s * (Beichmann et al. 1986| ). VLA observations show a weak, 
unresolved radio source with a flux density of ~ 6 mJy at 4.85 
GHz (Laurent-Muehleisen et al. 



1997) 



Wills et al. (1992) presented polarimetric and spectropho- 
tometric observations and found high polarization which rises 
from 1.4% in the K-band (2.2/i) to 8% in the U-band (0.36//). 
From the very large hydrogen broad-line ratios (H a /Hp = 
5.9; Pa a /H a = 0.15) and the shape of the optical continuum 
they derive E(B — V)~ 0.3, which corresponds to a hydrogen 
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1 Ho = 50 km s _1 Mpc -1 and qo = 0.5 are assumed throughout 
this paper. 



column density of 1.7 x 10 cm assuming standard dust and 
a Galactic gas-to-dust ratio. To explain their observations, Wills 



et al. (1992) argued for a model which invokes a bipolar geom- 



etry of the central region and in which the line-of-sight just 
grazes the edge of a dusty, parsec-scale torus or disk. The ob- 
served optical emission is a superposition of two components: 
a direct one, which gets attenuated as the light passes through 
the molecular torus, and a second one, which is emitted in polar 
directions and then scattered towards us by electrons or small 
dust grains. The latter component is therefore polarized parallel 
to the major axis of the host galaxy. 

IRAS 13349+2438 turned out to be an extraordinarily 
bright and variable soft X-ray source. In the ROSAT All-Sky 
Survey it was detected with m 2.5 cts/s in the 0.1-2.4 keV en- 



ergy band (Brinkmann & Siebert 1994). Brandt et al. (1996 



BFP96 hereafter) analyzed two short ROSAT PSPC pointed 
observations and found large amplitude variability as well as 
evidence for an ionized absorber in the soft X-ray spectrum 
by identifying the O VIII absorption edge at (rest frame) 0.871 
keV. They also suggested the presence of dust within the so- 
called 'warm absorber' in order to explain the apparent discrep- 
ancy between the observed optical reddening and the absence 
of any cold absorption in excess of the Galactic Nr value. 

Whereas the model of a warm absorber with internal 
dust has meanwhile been successfully applied to the X-ray 
spectra of several AGN (e.g. NGC 3227, NGC 3786, IRAS 



17020+4544; Komossa & Fink 1997a,b; Leighly et al. 1997b 



Komossa & Bade 1998) this has not yet been done for IRAS 
13 349+243 8 (for first results see Komossa & Greiner |1999[). 



Given the potentially strong modifications of the soft X-ray 
spectrum by the presence of dust (e.g. Komossa & Fink 1997a| ) 
it is important to scrutinize whether (and under which condi- 
tions) a dusty warm absorber is consistent with the observed 
X-ray spectrum. ROSAT spectra are best suited to accomplish 
this, since the most pronounced signature of dust is the carbon 
edge at 0.28 keV, which is outside the ASCA energy range. 



Brinkmann et al. (1996) analyzed the ASCA observation of 



IRAS 13349+2438 and found distinct changes in the spectral 
parameters compared to the earlier ROSAT data. These find- 



ings were confirmed by Brandt et al. (1997) in an independent 



analysis of the data. In particular, the strong O VIII absorption 
edge found in the PSPC observation was not detected anymore. 
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Fig. 1. Upper panel: the surface brightness profile (filled cir- 
cles) and the best-fitting PSF model for IRAS 13349+2438 
(solid line). Lower panel: ratio of surface brightness data to 
PSF model. 



Instead, both analyses find a spectral feature around ~ 0.65 
keV, which is interpreted as an emission line by Brinkmann et 
al. (1996 ) and a weak Ovil absorption edge by Brandt et al. 
dl997b . 

Besides a successful spectral fit, a key argument for invok- 
ing the presence of dusty warm absorbers is the apparent dis- 
crepancy between the column densities of cold material derived 
from the optical reddening and the absence of excess soft X-ray 
absorption. A crucial assumption for this argument is that both, 
the optical and the X-ray continua, are seen along the same 
line of sight. The necessity of dusty warm material would be 
alleviated if the X-ray emission were spatially extended, i.e. 
if the X-rays would not pass through the material causing the 
optical extinction. Observations with the ROSAT HRI allow to 
set limits on the contribution of an extended component to the 
X-ray emission by its improved spatial resolution and the con- 
firmation of the r apid X -ray variability observed with ASCA 
(Brinkmann et al. |1996|). 



We note that IRAS 13349+2438 was also considered as 



a narrow-line Seyfert 1 galaxy (Osterbrock & Pogge 1985 ), 
due to some of its optical (permitted emission line widths, 
[OlII]/H/3 ratio) and X-ray properties (steep spectrum, rapid 
variability), which are typical for this class of AGN. We will 
discuss some of our results in view of this hypothesis. 

2. The ROSAT HRI data 

IRAS 13349+2438 was observed twice with the ROSAT HRI. 
The first observation (hereafter HI) took place in 1996 from 
July 16 to July 21 and resulted in an effective exposure of 6209 



HI 



H2 



IB 17 18 19 20 21 22 15 16 17 IB 19 20 21 22 23 24 25 
Day of July 1996 Day of June 1997 

Fig. 2. The HRI light curve of IRAS 13349+2438 for both ob- 
servations. IRAS 13349+2438 displays large amplitude count 
rate variations on time scales of hours to days. The fastest vari- 
ability detected is a 20% decline in count rate within « 3 hours. 



s. The second observation (H2) was conducted roughly one 
year later, from June 15 to June 24 1997, and resulted in 24748 
s of good exposure. IRAS 13349+2438 was observed on-axis 
in both observations with an average count rate of w 0.3 cts/s. 

2.1. Spatial extent 

The high spatial resolution of the ROSAT HRI provides tight 
constraints on any extended emission component in IRAS 
13349+2438. The HRI point spread function (PSF) can be 
approximated by two Gaussians with <j\ w 2.2 arcsec and 
(72 ~ 4.0 arcsec and an exponential term, which is relevant 
for la rge radii (« 30 arcsec) and bright sources (David et al. 



1997[ ). The inner core of the PSF corresponds to a projected 
linear size of « 5.7 kpc at the redshift of IRAS 13349+2438. 
Due to random errors in the aspect solution, the values of <7i 
and 02 vary from 1.9 to 2.5 arcsec and from 3.3 to 4.1 arcsec, 
respectively, between individual observations of point sources. 
Further, elongations of the PSF by uncorrected residual wob- 
ble motion have been reported for a number of sources (Morse 



1994; W.Pietsch, priv. com.). 



In a practical approach, the observation dependent uncer- 
tainties were modeled by smearing the theoretical PSF with 
an additional Gaussian characterized by a+. We searched for 
the best-fitting a + by fitting the smeared PSF to the binned 
radial profile of IRAS 13349+2438 from the H2 observation 
using a least-squares method. The best description, shown in 
Fig. |lj is achieved for c+ ~ 1.8 arcsec. This value is well 
within the range determined from numerous other analyses of 
point source observations. Given the mentioned uncertainties 
in the determination of the correct HRI PSF for an individ- 
ual observation, we conclude that the X-ray emission of IRAS 
13349+2438 is point like to the limit of the HRI resolution. 
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Fig. 3. The historic X-ray light curve for IRAS 13349+2438. Shown are the All-Sky Survey observation (RASS), both PSPC 
pointings (P1,P2), the ASCA data and the two HRI observations (H1,H2). All measurements were converted to ROSAT PSPC 
count rates assuming a simple power law spectrum with Y = 2.8 and Galactic absorption. 



2.2. Light curves 

Source counts were extracted from a circular region with radius 
1.5 arcmin around the centroid position of the X-ray source. 
The background was determined from a source free annulus 
with inner radius 3 and outer radius 5 arcmin. Since IRAS 
13349+2438 was observed on-axis, no vignetting correction 
was applied to the data. Finally, the observations were split into 
individual observation intervals with durations between 1000 to 
3000 s. The intervals were chosen by hand to ensure sufficient 
photon statistics in each data point and to match as closely as 
possible the distribution of the observation intervals in "real 
time". The resulting light curves are shown in Fig. |[ 

The average count rate in both, HI and H2, is w 0.3 cts/s. 
Assuming a simple power law spectrum with a photon index 
of r = 2.8 (see Sect. 3.3.1) and Galactic absorption (JVh = 



(1.1 ± 0.2) x 10 20 cm -2 ; Murphy et al. 1996 1, this count rate 



yields an unabsorbed 0.1-2.4 keV flux of 8.2 x 10 12 erg cm 2 
s _1 . The corresponding rest frame 0.1-2.4 keV luminosity is 
4.6 x 10 44 erg s" 1 . 

IRAS 13349+2438 clearly displays rapid and large ampli- 
tude variability in both HRI observations. At the beginning of 
H2 we obviously observed the decline of a larger outburst and 
the count rate decreases by almost a factor of two within one 
day. The fastest variability seen is a decline by 20% within 
about 3 hours, which is truly remarkable given the high X-ray 
luminosi ty of IRAS 13349+2438. Following Lawrence & Pa- 
padakis (1993 1, who find a correlation between luminosity and 
doubling time scale from EXOSAT variability power spectrum 



analysis, a variability time scale of the order of 10 6 to 10 sec 
would have been expected for IRAS 13349+2438. 

For illustrative purposes we show in Fig. || the historical 
X-ray light curve for IRAS 13349+2438 including all available 
data from the ROSAT All-Sky Survey to the latest HRI obser- 
vations. All measurements have been converted to PSPC count 
rates assuming a simple power law spectrum (Y = 2.8) and 
Galactic absorption. 

X-ray variability was detected in all previous observations 
of IRAS 13349+2438. BFP96 note a factor of w 4 change be- 
tween the two PSPC observations separated by about one year 
(PI and P2). During P2 in December 1992, IRAS 13349+2438 
obviously happened to be in a historical high state. Unfortu- 
nately, P2 only lasted rj 1570 s and no variability was de- 
tected within this time interval (BFP96). In the ASCA observa- 
tion a 50% increase in flux within about 5 hours was detected 



(Brinkmann et al. |1996| ). This time scale is similar to that found 
in the HRI observation. 

Using the observed variability time scales and causality ar- 
guments we can derive upper limits on the size of the X-ray 
emitting region (R < c ■ At; neglecting relativistic beaming). 
A time scale of one day gives an upper limit of ks 180 AU. The 
fastest observed variability time scales of 3 to 5 hours suggest a 
size of the X-ray source not larger than 20 to 30 AU. This prac- 
tically rules out a significant contribution of scattered X-rays, 
since the observed sizes of electron scattering mirrors are a fac- 
tor of 10 4 - 10 5 larger (~ 50 pc in NGC 1068; Antonucci et al. 



1994J). In addition, a small mirror would imply very high elec- 



tron densities and a Thomson-thick scattering medium, which 



4 
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might be in contradiction to observations at other wavelengths 
(see discussion in BFP96). 

Variability also constrains the efficiency of the conversion 



of infalling matter into radiation energy (e.g. Fabian 1979) 

AL [erg/ s 



V 



> 



At[s 



2.1 x 10 42 

The observed change in luminosity of AL w 2.3 x 10 44 



erg 



s _1 within w 10 4 s gives 77 <; 1.2%. 

We finally note, that the above arguments hold only in 
the absence of relativistic beaming, which would effectively 
shorten the observed variability time scales by time dilation ef- 
fects. Models involving relativistic motion have recently been 
proposed for a number of narrow-line Seyfert 1 (NLS1) galax- 
ies to explain spectral features (Leighly et al. 1997a; Hayashida 



1997) and in particular to explain the rapid X-ray variability in 



IRAS 13224-3809 (Boiler et al. |1997[ ). Despite its high lumi- 
nosity, IRAS 13349+2438 shows many similarities to NLS1 
galaxies, e.g. a steep soft X-ray spectrum and rapid variabil- 
i ty. It al so satisfies the NLS 1 classification criteria of Goodrich 
( 1989b ). Therefore the presence of relativistic effects in IRAS 
13349+2438 cannot be excluded. 

3. The ROSAT PSPC data: Spectral analysis 

3.1. Data preparation 

IRAS 13349+2438 was observed in two short pointings with 
the ROSAT PSPC from January 7 to January 14, 1992 (t e ff « 
3280 s) and on December 16 1992 (t cS w 1570 s). Here- 
after the two observations are referred to as PI and P2, re- 
spectively. The data were extracted from the ROSAT archive at 
MPE and analyzed using standard routines within the EXSAS 
environment (Zimmermann et al. 1994). Source photons were 
extracted from a circular cell with radius 3 arcmin. The back- 
ground was determined from a circular source-free region near 
the target source. Corrections for vignetting and dead-time 
were applied. For the spectral analysis, the pulse height spec- 
trum between channels 1 1 and 240 was rebinned such that the 
signal-to-noise ratio in each bin was at least 10. 

In the ROSAT All-Sky Survey IRAS 13349+2438 was ob- 
served for about «365 s. We extracted source photons from a 
circle with radius 10 arcmin centered on IRAS 13349+2438. 
To ensure a similar exposure, the background was determined 
from a source free circle with the same radius, but displaced 
from IRAS 13349+2438 along the scanning direction of the 
satellite during the All-Sky Survey. Again the data were cor- 
rected for vignetting and dead time. In total, about 930 photons 
were accumulated for IRAS 1 3349+2438 during the Survey ob- 
servation. Finally, the pulse height spectrum between channels 
11 and 240 was rebinned according to a constant signal-to- 
noise ratio of 5.hl626.tex 



3.2. Warm absorber models 

To discuss the X-ray spectrum in terms of warm absorption we 
used warm absorber models calculated with the photoioniza- 
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Fig. 4. Two representative models of a dust-free (top panel) 
and a dusty (bottom panel) warm absorber to illustrate the in- 
fluence of dust. The range of the abscissa corresponds to the 
ROSAT energy band (0.1-2.4 keV). The respective model pa- 
rameters correspond to the best-fit values determined for the 
ROSAT PSPC observations. In the case of the dust-free warm 
absorber these are V = 2.2 and log N w = 22.5. The dusty 
warm absorber model is shown for log N w = log N opt = 21.2 
and r = 2.9. All model spectra have been corrected for Galac- 
tic cold absorption. The straight lines mark the respective in- 
trinsic power law spectra. See text for further details. 



tion code CLOUDY (Ferland 1993 ). For the general assump- 
tions of the models see Komossa & Fink (1997a,b). Briefly, the 
ionized material was assumed to be of constant density (log njr 
= 7) and of solar abundances (Grevesse & Anders 1989) in the 



dust-free case. A typical Seyfert IR to 7-ray spectral energy 
distribution (with a uv -x = —1-4 in the EUV) was assumed 
as ionizing continuum. In the warm absorber models including 
dust, the dust composition and grain size distribution were like 
those of the Galactic diffuse interstellar medium (Mathis et al. 



1977), and the chemical abundances were depleted correspond- 
ingly, if not mentioned otherwise. We focused on Galactic- 
ISM-like dust, because this was also assumed in the optical 
reddening estimates (e.g. Brandt et al. 1996). To test specific 
scenarios we also varied the dust properties, but these are not 
treated as additional free fit parameters. The warm absorber 
properties which are derived from X-ray spectral fits are the 
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ionization parameter U = Q / (4irr 2 niic) and the hydrogen col- 
umn density N w of the ionized material, where Q is the count 
rate of incident photons above the Lyman limit, r is the dis- 
tance between the central source and the warm absorber and 
riH is the hydrogen density. 

To illustrate the influence of dust on a warm absorber, we 
show in Fig. |] the model spectra for a dust-free and a dusty 
warm absorber model in the ROSAT PSPC energy range. In 
the dust-free case the observed spectrum, when it is described 
with a simple power law model, appears to be steeper than the 
intrinsic power law continuum because of a sequence of ab- 
sorption edges due to ionized oxygen and neon. If dust is added 
to the warm absorber, the observed spectrum changes dramat- 
ically. A carbon absorption edge appears at 0.28 keV due to 
the graphite species of dust and the absorption edges of ionized 
oxygen are less pronounced, because of the dust-induced lower 
ionization state of the warm absorber gas (Komossa & Bade 



1998). In particular, the O VIII edge is considerably weaker rel- 



ative to O VII (and almost completely disappears in the case of 
IRAS 13349+2438). Contrary to the dust-free warm absorber, 
the presence of dust leads to an apparent flattening of the X-ray 
spectrum in the ROSAT PSPC energy band. 

3.3. Results 

3.3.1. Pointed Observations 

First, we applied various 'standard models' to the data, e.g. a 
simple power law and a power law plus black body, to check 
whether the results of BFP96 could be reproduced. The residu- 
als of the power law fit to the PSPC data are displayed in the top 
panels of Fig. |. We obtain T = 2.80 ±0.15, N K = 1.10±g;^ 
1.55 for PI and T = 2.83 ± 0.11, N H = 



cm and xid 



(1.22 ± 0.20) x 10 2U cm" 2 and X 



2 

red 



1.72 for P2. Clearly, 



simple power law models are ruled out for both observations. 

In a second step, a dust-free warm absorber model was ap- 
plied. In contrast to the approach of BFP96 we use the addi- 
tional information on the hard X-ray power law available from 
the ASCA observation, i.e. we fix the intrinsic power law in- 
dex to Tasca 



1996 



see also Sect. 



2.2 (Brinkmann et al. 
4.2). We find an acceptable fit of the X-ray spectrum in terms 
of a dust-free warm absorber, with log N w = 22.5 ± 0.3 
(X 2 ed = 1.05) and log N w = 22.7 ± 0.2 (x r 2 cd = 0.84) for 
PI and P2, respectively (Fig. |^, middle panels). The value of 
JV W is consistent with that of BFP96. The amount of cold ab- 
sorption in this model is slightly below, but still consistent with 
the Galatic value (AT H = (0.08 ± 0.02) x 10 21 cm~ 2 , as com- 
pared to A^ H ,gai = 0.11 x 10 21 cm~ 2 ). 

We conclude that the dust-free warm absorber model pro- 
vides a successful fit to the ROSAT X-ray spectra and removes 
almost all systematic residuals. Moreover, the intrinsic X-ray 
continuum is consistent with the hard X-ray photon index ob- 
tained with ASCA. 

If we now modify this best fit model by adding the ap- 
propriate amount of dust to the warm absorber which corre- 
sponds to the value of N w , taking into account the gas-dust and 



gas-radiation interaction, the expected X-ray spectrum changes 
drastically and the data cannot be fit at all (x 2 od ~ 150). This is 
partly due to dust-induced additional absorption, in particular 
at the carbon edge, which leads to a much flatter spectrum than 
observed. This still holds if we allow for non-standard dust, i.e., 
if we selectively exclude either the graphite or silicate species. 
A more detailed exploration of the very large paramater space 
of dust properties is not warranted with current data and should 
await the availability of higher resolution X-ray spectra. 

The expected column density derived from optical extinc- 
tion is much smaller than the value of N w determined from the 
dust-free warm absorber fit. Therefore, in a next step, we al- 
lowed all parameters (except T = Tasca) to vary and checked 
whether a dusty warm absorber can be successfully fit at all. 
This is not the case. For example, if log N w is fixed to log 



N, 



opt 



21.2 we get X ; c 



40. 



Next, the underlying power law spectrum was allowed to 
vary as well. Fixing N w to A^pt again, the fit of the P2 observa- 
tion results in a steep intrinsic spectrum with T = 2.86 ± 0.04 
in order to compensate for the 'flattening effect' of dust (cf. 
Komossa & Fink 1997a,b). The intrinsic continuum is much 
steeper than obtained with ASCA and still the fit is statisti- 
cally not acceptable (x 2 od = 1.68). The ionization parameter 
is not well constrained by the data (log £7 = 0.19to'gg). Sim- 
ilar results are obtained for the PI observation (x 2 cd = 1-35). 
The residuals between 0.6 and 1 keV remain (Fig. g, bottom 
panels) and a dust-free warm absorber model provides a better 
description of the PSPC data. 

3.3.2. Survey observation 

A simple power law model provides an excellent fit to the Sur- 
vey data with T = -3.2 ± 0.3 and N H = (0.15 ± 0.06) x 10 21 
cm" 2 (x? cd = 0.74). 

When fitting the warm absorber models to the Survey ob- 
servation we get the following results: A dust-free warm ab- 
sorber with T = —2.2 provides a successful fit and the result- 
ing model parameters are consistent with those of the pointed 
observations within the errors. If we add dust to the model and 
fix Ny, to A^pt, no successful fit can be achieved. For example, 
fixing Ah to the Galactic value, the fit still yields an unaccept- 
ably high x 2 ed °f about 11. Allowing for a steeper intrinsic 
spectrum, with T = 2.9 as in the pointed observations, a dusty 
warm absorber with N w = Aopt is consistent with the spectrum 
from the Survey observation. However, the quality of the fit is 
slightly worse compared to the dust-free warm absorber. 

We conclude that, within the limits of low photon statistics, 
the Survey data are roughly consistent with the pointed obser- 
vations. 



4. The ASCA data 

The ASCA observation of IRAS 13349+2 438 w as analyzed 
and di scusse d in detail by Brinkmann et al. ( 1996 ) and Brandt 
et al. ( 1997 ). Therefore we describe in the following only the 
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Fig. 5. The residuals of a simple power law fit (top), a warm absorber (middle) and a dusty warm absorber fit to both PSPC 
observations. Left: PI, right: P2. 



differences to the above papers in terms of data preparation and 
then turn to our new spectral analysis. 

4.1. Data preparation 

We used the 'Revision 2' data obtained from the ASCA public 
archive at Goddard Space Flight Center (GSFC) and applied 
the following conservative screening criteria: For the GIS the 
minimum elevation angle above the Earth's limb (ELV) was 
chosen to be 5°. In the case of the SIS we used ELV> 10°. To 
avoid atmospheric contamination, data were only accepted in 
the SIS when the angle between the target and the bright earth 
(BR_EARTH) was greater than 20°. A minimum cut-off rigidity 
(COR) of 6 GeV/c was applied for both, the SIS and the GIS. 
Data taken within four read-out cycles of the CCD detectors 
after the passage of the South Atlantic Anomaly (S AA) and the 
day-night terminator are not considered in the analysis. Further, 
periods of high background were manually excluded from the 
data by checking the light curve of the observation. 

Source counts were extracted from a circular region cen- 
tered on the target with a radius of 6' for the GIS and 4' for 



the SIS. We used the local background determined from the 
observation in the analysis for both detectors. In particular, the 
GIS background was estimated from a source free region at the 
same off-axis angle as the source and with the same size as the 
source extraction region. 



Both Brinkmann et al. (1996) and Brandt et al. (1997) nei 



ther find spectral differences between the two SIS detectors nor 
between the two observations of IRAS 13349+2438 separated 
by about three days. We therefore merged the data from the 
two observations to increase the signal-to-noise ratio. The same 
procedure was applied to the GIS detectors. 

All spectra were rebinned to have at least 20 photons 
in each energy channel. This allows the use of the % 2 
technique to obtain the best fit values for the model spec- 
tra. We used the latest GIS redistribution matrices avail- 
able (V4_0) from the calibration database and created the 
SIS response matrices for our observation using SISRMG, 
which applies the latest charge transfer inefficiency (CTI) table 
(sisph2pi_l 103 97 .fits). The ancillary response files for 
all four detectors were generated using the ASCAARF program. 
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Spectra were fitted in the energy range 0.8 to 9 keV for the 
GIS and between 0.6 and 8 keV for the SIS. The upper en- 
ergy boundaries are given by the maximum energy at which 
the source was detected in each instrument. The lower energy 
boundaries result from the calibration uncertainties of the de- 



i SIS o GIS 2 
. SIS 1 a GIS 3 



Power low mode 



tectors (e.g. Dotani et al. 1997). If not mentioned otherwise 



we quote the results of simultaneous fits to all four instruments 
only. 

4.2. Spectral analysis 

To check for the consistency of our data preparation procedure, 
we first fitted several simple models to the ASCA dat a and 
compared them to previous analyses (Brinkmann et al 



1996 



Brandt et al. 1997[ Leighly 1999[ ). For the simple power law 
model we get 7Y H = (1.75±^ x 10 20 cm~ 2 , T = 2.22±0.05 
and x 2 = 548.3 (565 d.o.f). The fitted 7Y H is consistent with 
the Galactic value and is therefore fixed to 1.1 x 10 20 cm -2 
in the following. Restricting the fit to energies above 2 keV 



. SIS o GIS 2 
. SIS 1 a GIS 3 




Power law + edge mode 



we find T 



2.24 



+0.08 
-0.09 



and x 



187.8 (237 d.o.f.). In ac- 



cord with Brinkmann et al. (1996) and Leighly (1999), we note 
that the inclusion of a narrow {a = 0.05 keV) iron Ka line 
improves the fit significantly (A% 2 = —7.2; P = 99.9% ac- 
cording to an F-test for two additional parameters). The best- 



fitting parameters are: T = 2.28 ± 0.09, E Vmc = 6.8lt°;i 2 and 



EW = 3131 

204 The obtained iron line parameters are typ- 
ical for narrow-line Seyfert 1 galaxies (Leighly 1999). A gaus- 
sian emission line with the above parameters is included in all 
following models. Dividing the ASCA data by an extrapolation 
of this model to lower energies gives the residuals shown in the 
upper panel of Fig. |[ Clearly, there are systematic residuals 
below 1 .5 ke V. 

Adding a single absorption edge to the power law model 
results in the following best-fitting parameters: T = 2.33 ± 
0.04, E cdgc = 0.75 ± 0.04, r cdgc = 0.43_t£^ andx 2 = 509.5 




(563 d.o.f.). The agreement with the Brandt et al. ( |1997| ) results 
is excellent. As can be seen in the middle panel of Fig. the 
fit is very good and all systematic residuals have disappeared. 



As shown by Brandt et al. ( 1997 ), a dust free warm absorber 



Energy (keV) 

Fig. 6. The ratio of the ASCA data to various spectral models 
as described in the text. The four detectors are indicated by 
different symbols. The data have been rebinned according to 
the spectral resolution of the individual detectors. 



trum or whether calibration uncertainties of the SIS detectors 
are important, is difficult to answer. 

5. Discussion 



P0.07 



provides a good fit to the ASCA data with T = 2.31_ 15 
log N w = 21.59+°;^. 

Next we applied the dusty warm absorber model described 
in Sect. 3.2 to the ASCA data (including the iron Ka line with 
line energy and a fixed to the above values). The best-fitting 



and 5.7. A dust-free warm absorber? 



model parameters are as follows: T = 2.39_g :i, log U = 
0.00±£;^, log 7Y W = 21.33±g;|| and x 2 = 524.4 (563 d.o.f.). 
The ratio of data to model for this fit is shown in the lower 
panel of Fig. ||. Note that the column density of the ionized 
material corresponds very well to the value expected from the 
observed optical extinction. Although the fit is statistically ac- 
ceptable and clearly better than the simple power law model, it 
is significantly worse than the power law plus absorption edge 
fit. The data lie systematically above the model in the lowest 
energy channels. Whether this indicates that the dusty warm ab- 
sorber model is not an optimal description of the ASCA spec- 



As we have shown in Sect. 3, a dust-free warm absorber model 
represents the formally best description of the PSPC spectrum 
of IRAS 13349+2438. The systematic structures in the residu- 
als almost completely disappear and, in addition, the spectral 
shape of the intrinsic X-ray continuum is consistent with the 
hard X-ray spectrum as observed with ASCA (T « 2.3). We 
further note that applying the dust-free warm absorber model 
we also obtain a viable description of the ASCA data, although 
with significantly different physical parameters for the warm 
absorber as compared to the ROSAT PSPC results. The energy 
of the absorption edge is significantly lower in the ASCA data 
and consistent with O VII, whereas the PSPC data clearly favor 
O VIII. Furthermore, the optical depth of the absorption edge is 
much higher in the PSPC data and, consequently, the column 
density of the warm absorber derived from the ROSAT obser- 
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vations (log N w = 22.7 ± 0.2) is about an order of magnitude 
higher than in the ASCA observation (log N w = 21.59to26 ; 
Brandt et al. 1997). 

The rapid X-ray variability observed in IRAS 13349+2438 
indicates changes in the ionizing continuum and thus a vari- 
able ionization state of the warm absorber is not unlikely. This 
might explain the change in the absorption edge energy as 
IRAS 13349+2438 was about a factor of four brighter in X- 
rays during the second ROSAT observation (P2) as compared 
to the ASCA observation. In the bright state, most of the oxy- 
gen is ionized to O VIII. When the ionizing flux decreases, 
the O VIII ions recombine to O VII and the observed absorp- 
tion edge energy decreases. This kind of variation of the ion- 
ization state of the warm absorber has also been claimed for 
the Seyfert 1 galaxy MCG-06-30-15 (Reynolds et al 



1995 



Otani et al. 1996). The first PSPC observation (PI), when IRAS 



13349+2438 was in a state comparable to the ASCA observa- 
tion, suffers from low photon statistics and although the result- 
ing spectral parameters of the dust-free warm absorber model 
are consistent with P2 we note that the higher x 2 c( j (1-05 com- 
pared to 0.84) might indicate spectral changes. 

Apart from changes in the ionization state, the differences 
between the ASCA and ROSAT spectra also require a change 
in the column density of the ionized material of at least 4 x 10 22 
cm -2 between the two observations, i.e. within about 3 years. 
One might speculate that isolated clouds of ionized material 
stripped off the torus and moving across the line of sight could 
be responsible for such a change in column density. 

A principal difficulty of any dust-free warm absorber model 
of course is to explain the discrepancy between the observed 
amount of optical reddening and the absence of any cold X-ray 
absorption. Apart from postulating atypical gas-to-dust ratios 
in IRAS 13349+2438 we might think of two possible ways out: 
variable optical extinction and different paths for the optical 
and the X-ray radiation. 

With regard to the first possibility we note that the X-ray 
and the optical observations were not simultaneous and there- 
fore variable optical extinction cannot be excluded. In fact, 
changes in the reddening of emission lines in Seyfert 1.8 galax- 



ies on a time scale of years were reported by Goodrich (1989c , 



199C) and interpreted in terms of moving cold obscuring ma- 



terial. On the other hand, no changes in the optical extinction 



have been observed up to now (Wills et al. 1992 ; Lanzetta et al 



1 993 ) and the preliminary results of recent spectroscopic obser- 



vations of IRAS 13349+2438 (Leighly, priv. com.; Papadakis, 
pri v. com .) are completely consistent with the original Wills et 
al. ( 1992 ) results. Furthermore, the amount of cold absorption 
measured with the ROSAT PSPC and ASCA is always con- 
sistent with the Galactic value and does not indicate variable 
cold absorption. Both results favor constant optical extinction. 
Obviously, simultaneous soft X-ray and optical observations of 
IRAS 13349+2438 are required to ultimately decide whether 
variable optical extinction plays a role, but based on current 
observations this seems highly unlikely. 

An alternative explanation for the apparent discrepancy be- 
tween optical and cold X-ray absorption would be a special 



geometry, in which the optical and the X-ray continuum do not 
travel along the same paths to the observer. Since the observed 
rapid X-ray variability of IRAS 13349+2438 effectively ex- 
cludes a significant contribution of scattered X-rays, one might 
invoke a partial covering geometry, where the optical contin- 
uum happens to pass through a lower total column density of 
obscuring material. Although such a scenario might be con- 
structed for a single object, this approach becomes increasingly 
contrived in the light of the growing number objects, for which 
a 'simple' dusty warm absorber provides a consistent explana- 
tion of the X-ray and optical properties (e.g., NGC 3227, NGC 
3786, IRAS 17020+4544). We therefore now turn to the dusty 
warm absorber hypothesis. 

5.2. A dusty warm absorber? 

The prime argument for the presence of dust within the warm 
absorber is the apparent discrepancy between the observed op- 
tical extinction and the absence of cold absorption in the X- 
ray spectrum. Furthermore, as we have shown in Sect.4.2, a 
self-consistent Galactic-ISM dusty warm absorber model is 
a reasonable representation of the ASCA spectrum of IRAS 
13349+2438 and the column density of the warm absorber de- 
termined by spectral fitting is consistent with the A^h derived 
from optical reddening. 

We note, however, that the ASCA spectrum is not very sen- 
sitive to the dustiness of the warm absorber, since the most im- 
portant features (like the Carbon absorption edge) are expected 
below 0.6 keV. Nevertheless, the absence of a strong Ovm 
edge in the ASCA spectrum^ (Brinkmann et al. 1996) is ex- 
actly the expected signature of dust within a warm absorber 
(cf. Fig. Q). The structure in the spectrum around 0.65 ke V was 
modeled as an Ovil absorption edge by Brandt et al. ( 1997| ) 
and as an isolated emission line by Brinkmann et al. ( 1996 ). 
In fact, even both features may be present and overlapping. We 
note that also an instrumental emission feature has been re- 



ported around 0.6 keV (see Sect. 2. 1 of Otani et al. |1996| ). 

The major problem with dusty warm absorber models for 
IRAS 13349+2438 lies in the fact that they are, at most, 
marginally consistent with the ROSAT PSPC data. Even if one 
accepts the reduced quality of the fit (x 2 d ~ 1-7), a temporal 
change of the intrinsic continuum of Ar ks 0.5 — 0.7 between 
the ROSAT and the ASCA observation has to be explained. 

The difference in the quality of the fit between the dust-free 
and the dusty warm absorbers is due to the clear indication of 
an O VIII absorption edge in the ROSAT data (see also Table 1 
of BFP96), which should be much weaker if the warm gas is 
dusty. In principle, it is possible to introduce increased O VIII 
absorption also in the case of a dusty warm absorber, for exam- 
ple by increasing the gas-phase oxygen abundance (which is 
depleted in the models to account for the binding of metals into 
dust). However, the limited spectral resolution of ROSAT cer- 
tainly does not warrant this kind of fine-tuning to achieve opti- 



2 This is confirmed by Brandt et al. (1997 ), who find r = Omtlm 



for the O VIII absorption edge, which is consistent with zero. 
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mal fits and the improved spectral capabilities of future X-ray 
missions are be needed to study in detail the elemental abun- 
dances of the warm absorber as well as the properties of the 
any dust in IRAS 13349+2438. 

As far as the change in the intrinsic continuum is con- 
cerned, one might speculate about various explanations, includ- 
ing a 'real' change in the spectrum and instrumental effects. 

The X-ray continuum spectrum might indeed have varied 
by AT w 0.5 - 0.7 between the ROSAT and the ASCA obser- 
vations, which are separated by about three years. This possi- 
bility of course cannot be ruled out and in fact spectral changes 
have previously been observed in narrow-line Seyfert 1 galax- 
ies (e.g. Mrk 766; Leighly et al. 1996). However, we consider it 
unlikely for IRAS 13349+2438 in view of the fact that no spec- 
tral variations have been observed between as well as within 
the individual PSPC and ASCA observations. In particular, the 
spectra observed in the two PSPC observations seem to be con- 
sistent despite a dramatic increase in flux by a factor of four. 
Furthermore, no spectral variations have been detected within 
the ASCA observation (e.g. Brinkmann et al. 1996). 



Systematic differences between ASCA and ROSAT spectra 
have been reported for various simultan eous observations in 
the past (e.g. NGC 5548; Iwasawa et al. ( 1999 ) and references 
therein). In particular, the photon index tends to be steeper by 
AT w 0.4 in these PSPC observations as compared to ASCA. 
On the other hand there are also observations, where the ASCA 
and ROSAT spectra do agree reasonably well (e.g. Miyaji et 
al. |1998[ Brinkmann et al. |1998[ Cappi et al. |1997|). Although 



these observations were not done simultaneously, it seems un- 
likely that spectral changes occured in all these sources which 
exactly compensate for any putative calibration uncertainties. 
The ASCA/ROSAT discrepancies definitely need to be further 
investigated. This is, however, beyond the scope of this paper. 

We conclude that it cannot be excluded that the observed 
difference between the ROSAT and ASCA power law contin- 
uum spectrum of IRAS 13349+2438 is at least partly due to in- 
strumental effects. We note, however, that even if we allow for 
the steeper PSPC spectrum in the dusty warm absorber model 
(and only in this model a systematically steeper continuum is 
required for the ROSAT data), the fit is still only marginally 
acceptable. 

In the previous section we discussed the possibility of a 
variable dust-free warm absorber. In order to preserve the dusty 
warm absorber hypothesis, and hence to obtain a physical 
model to explain all X-ray and the optical properties consis- 
tently, we might also speculate about more complicated warm 
absorber models. For example, a two-zone warm absorber has 
been proposed for the Seyfert 1 galaxy MCG-6-30-15 (Otani 
et al. 1996[ ), i.e. a variable inner warm absorber responsible 
for the O VIII absorption edge and an outer warm absorber in 
a lower ionization state and thus mainly imprinting the O VII 
edge. IRAS 13349+2438 might be a similar case: the inner 
warm absorber varied between the ROSAT and the ASCA ob- 
servation either by changes in the ionization state and/or by 
a variation in N w , whereas the outer dusty warm absorber is 
present in both observations. 



6. Conclusions 

We presented new ROSAT HRI data for IRAS 1 3349+2438 and 
a re-analysis of the ROSAT PSPC and ASCA spectra in terms 
of a warm absorber, but self-consistently including the effects 
of dust on the soft X-ray spectrum. 

The HRI light curve of IRAS 13349+2438 confirms the 
rapid and large amplitude v ariabil ity already noted in the 
ASCA data (Brinkmann et al. 1996 1. We observed a factor of 



two change in count rate within one day and 20% variability 
within about three hours. These rapid variations of the soft X- 
ray flux rule out a significant contribution of scattered X-rays 
to the total flux. 

The analysis of the measured ROSAT PSPC and ASCA 
spectra gives contradictory results. The PSPC data are best fit 
by a dust-free warm absorber, but large changes in the column 
density of the ionized material and the ionization state are re- 
quired to describe the ASCA data with this model. And finally 
it cannot explain the discrepancy between the optical reddening 
and the absence of cold X-ray absorption. On the other hand, 
a dusty warm absorber is consistent with the ASCA data and 
gives a column density of the warm material, which agrees very 
well with the one predicted from optical reddening. However, 
the dusty warm absorber is, at most, marginally consistent with 
the PSPC data and requires a much steeper (AT w 0.5 — 0.7) 
intrinsic continuum as that found in the ASCA data. 

Physical effects, such as a variable warm absorber column 
density and ionization state as well as instrumental effects have 
to be considered to explain the apparent discrepancies. Maybe 
a two-zone warm absorber, similar to the one postulated for the 
Seyfert 1 galaxy MCG-6-30-15, is able to provide a consistent 
explanation of all X-ray data and the optical extinction, but we 
consider it an over-parameterization of presently available X- 
ray data. From spectral fits to current ROSAT and ASCA data 
alone the presence of dust within the warm absorber cannot 
be proven and the discrepancy between the observed optical 
reddening and the absence of cold X-ray absorption remains 
the strongest argument in favor of a dusty warm absorber in 
IRAS 13349+2438. Future X-ray observations (with XMM or 
probably SAX), will be able to directly probe the dust content 
of the warm absorber, e.g. via the predicted Carbon absorption 
edge. 
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